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ABSTRACT 

Context. We present an XMM-Newton analysis of the X-ray spectra of 39 clusters of galaxies at 0.4 < z < 1.4, covering a temperature 
range of 1.5 < kT < llkeV. 

Aims. The main goal of this paper is to study how the abundance evolves with redshift not only by means of a single emission measure 
performed on the whole cluster but also by spatially resolving the cluster emission. 

Methods. We performed a spatially resolved spectral analysis, using Cash statistics and modeling the XMM-Newton background 
instead of subtracting it, by analyzing the contribution of the core emission to the observed metallicity. 

Results. We do not observe a statistically significant (> 2<x) abundance evolution with redshift. The most significant deviation from 
no evolution (at a 90% confidence level) is observed by considering the emission from the whole cluster (r < 0. 67500), which can 
be parametrized as Z oc (1 + z )-° 8 ± 5 . Dividing the emission into three radial bins, no significant evidence of abundance evolution 
is observed when fitting the data with a power law. We find close agreement with measurements presented in previous studies. 
Computing the error-weighted mean of the spatially resolved abundances into three redshift bins, we find that it is consistent with 
being constant with redshift. Although the large error bars in the measurement of the weighted-mean abundance prevent us from 
claiming any statistically significant spatially resolved evolution, the trend with z in the 0. 15-0.4r 50 o radial bin complements nicely 
the measures of Maughan et al., and broadly agrees with theoretical predictions. We also find that the data points derived from the 
spatially resolved analysis are well-fitted by the relation Z(r,z) = Z (l + (r/0. 15r 50 o) 2 )~ a ((l + z)/1.6)~ y , where Z = 0.36 ± 0.03, 
a = 0.32 ± 0.07, and y - 0.25 ± 0.57, which represents a significant negative trend of Z with radius and no significant evolution with 
redshift. 

Conclusions. We present the first attempt to determine the evolution of abundance at different positions in the clusters and with 
redshift. However, the sample size and the low-quality data statistics associated with most of the clusters studied prevents us from 
drawing any statistically significant conclusion about the different evolutionary path that the different regions of the clusters may have 
traversed. 

Key words. Galaxies: clusters: intracluster medium - X-rays: galaxies: clusters 

1. Introduction mainly via thermal bremsstrahlung and is easily detectable in 

the X-rays. Although the amount of energy supplied to the intra- 

g , Clusters of galaxies are the largest virialized structures in the cluster me dium (ICM) by SNe explosions depends on several 

Universe, arising from the gravitational collapse of high peaks factors (e g the physical condition of the ICM at the epoch of 

of primordial density perturbations. They represent unique sign- the enrichment) and cannot be obtained directly from X-ray ob- 

posts where the physical properties of the cosmic diffuse baryons ser vations, the radial distribution of metals in the ICM, as well 

can be studied in great detail and used to trace the past history of as their abundance as a function of time, represent the 'footprint' 

cosmic structure formation (e.g.lPeebles 1993; Coles & Lucchin| of cosmic star formation history and are crucial to trace the ef- 

119951 |Peacock||l999|; |Rosati et al.||2002t |Voit||2005|). The grav- fect of SN feedback on the ICM at different cosmic epochs(e.g. 

itational potential well of clusters is permeated by a hot thin |Ettorill2005l : iBorgani et al.ll2008h . 

gas, which is typically enriched with metals ejected form su- Seyeral studies of ' ±e radial distributions of metals 

pernovae (SNe) explosions through subsequent episodes of star fa ^ i CM at low red shift are present in the litem- 

formation (e.g. lMatteucci & Vetto lam 1988; |Renzini|1997|) . This mre (e g Finoguen o7CTl20^a lDe"G randi & Molendi 2001 



gas reaches temperatures of several 10 7 K and therefore emits |f rwin & B regman 2001; Tamura et al] l2004t IVikhlininSl 

120051: iBaldi et alJ l2007t iLeccardi & Molendil l2008al) . On the 

Send offprint requests to: A. Baldi, e-mail: other hand just a handful of works have tackled the evolu- 

alessandro . balcli@oabo . inaf . it tion of metal abundance in the ICM at larger look-back times. 
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Measurements of the metal content of the ICM at high-z has 
been obtained with single emission-weighted estimates from 
Chandra and XM M-Newton ex posures of 56 clusters at 0.3 $ 
Z £ 1.3 in lBalestraet~ai](l2007l) (hereafter BLS07). They mea- 
sured the iron abundance within (0. 15-0.3)7?,,,,- and found a neg- 
ative evolution of Zp e with redshift, for clusters at z ^ 0.5 that 
have a constant average Fe abundance of « O.25Z , while objects 
in the redshift range 0.3 5; z ^ 0.5 have a Zf e tha t is significantly 
higher («; 0.4Z o ). This result was confirmed by Maughan et al. 
(2008) (hereafter MAU08) for a sample of 1 16 Chandra clusters 
at 0.1 < z < 1.3. Stacking the cluster spectra extracted within 
7-500 m different redshift bins, they found that the metal abun- 
dance in the ICM drops by a factor of ~ 50% between clusters 
in the local Universe and clusters at z — 1. This evolution is 
still present when the cluster cores (r < 0.15r5oo) are excluded 
from the abundance measurements, indicating that the evolution 
is not simply driven by the disappearance of relaxed, cool-core 
clusters (which are known to have enha nced core metal abun - 
dances) from the popu lation at z > 0.5 ( VikMi rdn et al J 120071) . 
lAnderson et alj d2009l) found a similar drop in abundance be- 
tween z = 0.1 and z ~ 1, by adding the data from the MAU08 
samp le and from a low r edshift (z < 0.3) XMM-Newton sam- 
ple dSnowden et al.ll2008l) . to an additional sample of 29 galaxy 
clusters at 0.3 < z < 1.3 observed by XMM-Newton. 

In this work, we aim to study the evolution of metal abun- 
dance with redshift by also considering its radial dependence 
(using two or three spatial bins) not only to disentangle the evo- 
lution of the cores from the rest of the clusters but also to study 
the evolution in the outer regions (r > 0.4r5oo). The plan of 
the paper is the following. In § 2, we describe the sample and 
the data reduction procedure. We describe our spectral analysis 
strategy in § 3 with particular focus on the background treatment 
procedure (§ 3.3). In § 4, we present the results obtained in our 
analysis mainly regarding abundance evolution with redshift. We 
discuss our results in § 5 and summarize our conclusions in § 6. 

We adopt a cosmological model with Hq = 70 km/s/Mpc, 
Sl m = 0.3, and = 0.7 throughout. Confidence intervals are 
quoted at lcr unless otherwise stated. 

2. Sample selection and data reduction 

In Table Q] we present the list of galaxy clusters analyzed in this 
paper. The selected clusters are all the clusters at z > 0.4 ob- 
served at least once by XMM-Newton and with sufficient signal- 
to-noise ratio (S/N) to allow a spatially resolved spectral analysis 
(at least in two spatial bins). To be included in the analysis, each 
bin was required to have a S/N larger than 60% and to have at 
least 300 MOS net counts. This was a trade-off in order to have 
both a sizeable sample and to avoid the introduction of excessive 
noise in the measure of the abundance. Figure Q] shows the red- 
shift distribution of the sample, which is clearly peaked at low 
redshifts with more than half of the sample being at z < 0.6. 

The observation data files (ODF) were processed to produce 
calibrated event files using the XMM-Newton Science Analysis 
System (SAS v9 . . 0) processing task EPPROC and EMPROC 
for the pn and MOS, respectively. Unwanted hot, dead, or flick- 
ering pixels were removed as were events due to electronic noise. 
However, the pn data were only used in the imaging analysis, 
because of the possible effects of an under- or over-estimation 
of the pn particl e background on the analys is of the extended 
sources (see e.g. iLeccardi & Molendill2008bl and XMM-ES AS 
cookboolfl). As we also tested on a representative subsample of 

1 ftp://xmm.esac.esa.int/pub/xmm-esas/xmm-esas.pdf 



our clusters (in redshift, flux, and temperature), this leads to in- 
consistencies with the MOS detectors in the measurements of 
temperature and abundance. On the basis of the systematic dif- 
ferences between MOS and pn, a joint fit of pn spectra with MOS 
spectra would be unjustified, since it would introduce both a bias 
in the best-fit values and artificially smaller errors (caused by the 
worsened spectral fits). 

To search for periods of high background flaring, light curves 
for pattern=0 events and for pattern < 4 in the 10-12 keV are 
produced for the pn and the two MOS, respectively. The soft pro- 
ton cleaning was performed using a double-filtering process. We 
extracted a light curve in 100s bins in the 10-12 keV energy band 
by excluding the central CCD, applied a threshold of 0.20 cts s , 
produced a GTI file, and generated the filtered event file accord- 
ingly. This first step allows most flares to be eliminated, although 
softer flares may exist such that their contribution above 10 keV 
is negligible. We then extracted a light curve in the 2-5 keV band 
and fit the histogram obtained from this curve with a Gaussian 
distribution. Since most flares had been rejected in the previous 
step, the fit was usually very good. We calculated both the mean 
count rate, fi, and the standard deviation, cr, applied a threshold 
of fi + 3cr to the distribution, and generated the filtered event 
file. The intervals of very high background were removed using 
a 3cr clipping algorithm and the light-curves were then visually 
inspected to remove the background flaring periods not detected 
by the algorithm. In Table [1] we list the resulting clean exposure 
times for the pn and the MOS detectors. 

We selected events with patterns to 4 (single and double) 
for the pn and with patterns to 12 (single, double, and quadru- 
ple) for the MOS. Further filtering was applied to remove events 
with low spectral quality. For the pn in particular, we removed 
the out-of-field events ((FLAG & 0x10000) = 0)), and the events 
close to bad pixels((FLAG & 0x20) = 0), on offset columns 
((FLAG & 0x8) = 0), and close to a CCD window ((FLAG & 
0x4) = 0). We merged the event files from each detector (and in 
some cases from multiple observations) to create a single event 
file. The total event file was used to extract a 0.5-8 keV band 
image of each cluster, but not for the spectral analysis where the 
spectra were extracted from each individual detector (and obser- 
vation) and fitted simultaneously (See §|3]below). The use of the 
pn data in the resulting image was motivated by the possibility of 
detecting the point sources to remove from the cluster spectra at 
a lower flux limit. A list of point-like sources was created using 
the EBOXDETECT task and then visually inspected to check for 
false detections. When one or more point sources were located 
inside an extraction region, circular regions of radius ^ 15" (de- 
pending on the source brightness), centered at the position of 
each point source, were excised from the spectral extraction. 

3. Spectral analysis strategy 

Our strategy for performing a spatially resolved spectral analysis 
involves the determination of the overdensity radius rsoo- This 
is necessary because we aim to study the cluster abundances 
and temperatures in annuli with inner and outer radii depending 
on the physical properties of the cluster instead of the number 
of counts in each annulus. We analyze each cluster into two 
spatial bins: the region of the cluster core (corresponding to 
< r < 0.15r5oo) and the region immediately surrounding the 
core (0.15r5oo < r < 0.4r5oo). For most clusters, a third bin 
at r > 0.4r5oo, considering the outskirts of the cluster, could 
also be analyzed. The spatial extension of the third bin was 
determined by considering concentric annuli of thickness 0. 1 rsoo 
at increasing radii (0.4r 50 o < r < 0.5r 5 oo, 0.5r5oo < r < 0.6r 50 o, 
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Table 1. XMM-Newton archival observations of distant clusters of galaxies. 
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\sl L4.3 Ll^l 




Obs. Date 


Obs. ID 


T 

1 cleanMOS 1 


T 

1 clean MOS2 


1 clean, pn 




Rjidicil bins' 7 










(ksec) 


(ksec) 


(ksec) 


(10 20 cm" 2 ) 
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U.4U / 


oaaa a a 
ZUUU INOV uo 


a 1 aa/1 An 1 n 1 
U1U040U1U1 


A 1 1 
41. / 


41.1 


^A T 

3U.3 


1 n 
1 ,u 




KAJ LZLj.j+VjZjj 


A 1 AO 


OAA1 rv>, . T A 

zOUl Dec 30 


nnc i 1 aq a 1 
005l34Uo01 


Zl.o 


OO A 


1 A 1 

14.3 


1 o 
1 .5 


z 


KXLJUojd.I+j /jo 


A A 1 1 
U.41 t 


OAAC Al . . i 1 A 

ZUU3 Uct 1U 


nim^c i qa 1 
030z38t5Ul 


O/l o 

Z4.z 


11 A 

z3.y 


1 A 1 

14. / 


^2 O 

5.1 


1,Z,J 


RXJ2228. 6+2037 


0.412 


2003 Nov 18 


0147890101 


24.6 


24.4 


19.3 


4.3 


1,2,3 


nvni aai a i /i 
KXCJ 1003.0+3234 


A /I 1 A 

0.410 


OAAC XT™- ni 

2003 NOV 03 


A " 1 A"> C O 1 AA 1 

0302381601 


OC A 

23.0 


OC 1 

23.1 


1 C "7 

13. / 


1 . / 


1 O Q 

1,2,3 


MS0302. 3+1 / 1 / 


A A 

0.425 


OAAO A >.™ TI 

2002 Aug 23 


m ni a a i a i 
01 12190101 


1 o c 
12.3 


1 o 1 

12.1 


1.5 


A A 

9.4 


1 o 

1,2 


KXCJ 1200.2-0848 


A /I /I A 

0.440 


OAAO r\,. , , AA 

200 / Dec 09 


Arm i lf\A A1 

0302430401 


OA 1 

2y.i 


OO o 

28.8 


OA A 

20.6 


5. 1 


1 O 1 

1,2,3 


IRAS09 104+4 109 


0.442 


2003 Apr 27 


0147671001 


12.1 


12.3 


8.3 


1.4 


1,2,3 


RXJ0221. 1+1958 


0.450 


2005 Jul 14 


0302581301 


5.4 


6.5 


1.2 


9.1 


1,2 


RXJ1347.5-1145 


0.451 


2002 Jul 3 1 


0112960101 


32.6 


32.4 


27.5 


4.9 


1,2,3 


RXJ131 1.5-0551 


0.461 


2006 Jan 17 


0302582201 


25.1 


26.0 


19.1 


2.4 


1,2,3 






2007 Jul 26 


at ai ioni ni 

0502430101 


34.6 


36.0 


21.2 






RXJ0522.2-3625 


0.472 


2005 Aug 14 


0302580901 


19.2 


19.4 


16.0 


3.6 


2 


RXJ2359. 5-3211 


0.478 


2005 Jun 13 


0302580501 


38.8 


38.5 


29.5 


1.2 


1,2 


CLJ0030+2618 


0.500 


2005 Jul 06 


0302581101 


15.0 


14.0 


- 


3.7 


1,2,3 






2006 Jul 27 


0402750201 


27.0 


27.4 


19.6 










2006 Dec 19 


0402750601 


28.5 


28.6 


20.6 






W ARP J 2 1 46 . 0+0423 


0.531 


2005 Jun 10 


0302580701 


20.7 


20.8 


17.7 


4.8 


2 


RXJ0018.8+1602 


0.541 


2007 Dec 14 


0502860101 


41.1 


42.7 


30.1 


3.8 


1,2 


MS0015. 9+1609 


0.541 


2000 Dec 29 


0111000101 


30.6 


30.1 


20.4 


4.0 


1,2,3 






2000 Dec 30 


ni 1 i nnm n i 

01 1 1000201 


5.5 


5.4 








WARPJ0848. 8+4456 


0.543 


i n n i /~\ tic 

2001 Oct 15 


aao c 1 cn i n i 

0085150101 


39.1 


39.6 


33.7 


2.8 


2 






ZUU1 Uct Zl 


aao; 1 CA^AI 


23.3 


25.2 


13.0 










2001 Oct 21 


0085150301 


25.2 


26.4 


17.9 






A 1 ! t i i c a aoo i 
CLJ 1354-0221 


0.546 


OAAO T. 1 A 

2002 Jul 19 


Ul IzzjUIUI 


6.8 


7.7 


0.8 


3.2 


2 


\i / a Dim i in a i a/' ~) ii n 
WARPJ 1419.9+0634L 


0.549 


OAAC T.,1 A"7 

2003 Jul 0/ 


AOA1 / 1(\ 1 n 1 

0303670101 


43.1 


43.4 


29.9 


2.2 


1,2 


A/Tca/ici a aqac 
MS043 1.6-0303 


A CCA 

0.330 


OAA/1 C~.~ 1 A 

2004 Sep lo 


AO AC A "7 A 1 A 1 

02056/0101 


oc c 

23.3 


OA 1 

2o.l 


1 A A 

iy.4 


Q A 

3.9 


1 O 1 

1,2,3 


RXJ0018. 3+1618 


0.551 


2000 Dec 29 


0111000101 


30.5 


29.6 


22.4 


3.9 


2 


WARPJ 1419.3+0638 


0.574 


r\ f\f\r 1 1 f\—i 

2005 Jul 07 


nom/"'7ni ni 

0303670101 


43.1 


43.4 


29.9 


2.2 


1,2,3 


MS2053. 7-0449 


0.583 


r\ f\f\ 1 XT 11 

2001 Nov 14 


n i 1 ^ 1 c\f\/~ n i 

01 12190601 


16.4 


16.5 


8.9 


4.6 


2 


MACSJ0647.7+7015 


0.591 


2008 Oct 09 


0551850401 


52.9 


53.7 


32.3 


5.4 


1,2,3 






<~\r\r\(\ w x f\ a 

2009 Mar 04 


ncci o c 1 mi 

0551851301 


33.2 


34.6 


18.4 






/"IT I | | -| A , I'll O 

CLJ ll 20+43 18 


0.600 


2001 May 08 


n i m o /" f\f\ n i 

0107860201 


18.3 


18.7 


13.6 


3.0 


1,2,3 


CLJ 1334+5031 


0.620 


2001 Jun 07 


ni i 1 i aai n 1 

Oil l 160101 


36.1 


36.5 


27.0 


1.1 


1,2,3 


MACSJ0744. 9+3927 


0.698 


2008 Oct 17 


Am otai ni 

0551850101 


40.4 


41.4 


22.4 


5.7 


1,2,3 






OAAA AJfn~ O 1 

2009 Mai 21 


AC CIOCIOAI 

0331531201 


/Q O 

63.8 


AT A 
6/.0 


Q C O 

35.8 






WARPJ 1 342. 8+4028 


0.699 


2002 Jun 08 


0070340701 


33.3 


33.6 


25.0 


0.8 


2,3 


CLJ 1103.6+3555 


0.780 


200 1 May 15 


0070340301 


7.1 


7.4 


4.4 


2.6 


2 






OAA/1 TVfl,.., 1 1 

2004 May 1 1 


AO ACQ-7A 1 A 1 

02033 /0101 


1 c o 
33.2 


QA A 
36.0 


OT A 






A/IC1 1 17 C i A A C 

Mo 1 13 /. 5+6623 


0. /82 


OAAA / 1 , , t AC 

2000 Oct 03 


AAA/1 O AAO A 1 

0094500201 


1 o o 
12.5 


1 A O 

lo. 2 


A A 

4.y 


1 A 
1 .0 


1 O 1 

1,2,3 


CLJ1216. 8-1201 


0.794 


2003 Jul 06 


0143210801 


24.0 


24.2 


19.7 


3.3 


2,3 


^ /tf 1 AC ,1 /I ATM 

MS 1054.4-0321 


0.823 


2001 Jun 21 


AAA A O AA 1 A 1 

0094800101 


22.7 


23.6 


19.1 


3.1 


1,2,3 


CLJ0132. /-133 / 


A O 1 1 

0.831 


oaao r\,. O/l 
2002 Dec 24 


A 1 AAC A A 1 A 1 

0109340101 


CA A 

30.6 


CA A 

30.9 


Art i 
40.3 


1.3 


1 O Q 

1,2,3 


CLJ1226.9+3332 


0.890 


2001 Jun 18 


0070340501 


10.3 


10.7 


5.3 


1.8 


1,2,3 






2004 Jun 02 


0200340101 


67.2 


67.5 


53.4 






CLJ 1429.0+4241 


0.920 


2006 Jan 08 


0300140101 


33.1 


33.7 


18.3 


1.1 


1,2 


XLSSC029 


1.050 


2005 Jan 01 


0210490101 


82.0 


83.2 


64.3 


2.3 


2,3 


RDCSJ1252-2927 


1.237 


2003 Jan 03 


0057740301 


66.1 


66.2 


54.9 


6.1 


2,3 






2003 Jan 11 


0057740401 


65.4 


66.7 


56.1 






1WGAJ2235. 3-2557 


1.393 


2006 May 03 


0311190101 


76.0 


76.9 


61.1 


1.5 


2 



Notes. ia) Radial bins considered in the spatially resolved spectral analysis for each cluster: 1 = (0 - 0. 15r5oo), 2 = (0. 15 - 0.4rsoo), 3 = (> 0.4r5oo). 

etc.). We computed the source and background counts in each 
annulus moving towards the cluster outskirts, and we ceased 
to add annuli to the third bin spectrum when the source counts 
in the annulus fell below 30% of the total counts. We decided 
to include in the analysis only the spectra with at least 300 
net counts, in order to avoid the introduction of excessive 
noise in the measure of the abundance. In Table [TJ the radial 
bins considered in the spatially resolved spectral analysis 
are indicated. We have a total of 26, 39, and 24 spectra in the 
0-0.15r5oo, 0.15 — 0.4rgoo, and > 0.4rsoo radial bin, respectively. 



3.1. Determination ofr 5Q0 

As an implication to the M - T relation, the overdensity radius 
rsoo scales with the temperature of the cluster. To co mpute rsoo, 
we adopted the formula derived by Vikhlinin] (I2006I) 

r 500 /i£(z) = 0.792 Mpc, (1) 

where E(z) = (Q m (l + z) 3 + ^a) 5 . 
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Table 2. X-ray global properties. 



Cluster 


<r) 


(Z) 


r 500 


To-0.6r 500 


Zo-0.6r 50 o 




(keV) 


(Z ) 


(Mpc) 


(keV) 


(Z G ) 



A851 5.66 + 0.27 0.26+^ 0.975 5.73!°|! 0.24 + 0.06 

RXJ1213.5+0253 < 0.16 0.532 2.90!™ - 25 -om 

RXCJ0856. 1+3756 6.23 ± 0.57 0.42 + 0.11 1.023 6.43^ 0.37 ± ot09 

RXJ2228.6+2037 7.60 + 0.42 0.28 + 0.06 1.137 7 -89!° : | 0.28 ± 0.05 

RXCJ1003.0+3254 2.96^;|| 0.24^-^ 0.689 3 - 27 !nia 0.52 ±0.12 

MS0302.5+1717 3.97+° : " 0.04^^ 0.799 4.36™ : || 0.17+°J| 

RXCJ1206.2-0848 H-Ol^S 0.26 + 0.06 1.361 10.48! ,* 1 0.27 + 0.04 

IRAS09104+4109 4.91+°| 0.4l+° \ 3 2 0.886 5.49!";' 0.53+0.08 

RXJ0221. 1+1958 6.97!' : jj? < 0.29 1.062 5.96!™ 0.01!°^ 

RXJ1347.5-1145 11.23 ±0.43 0.23 ± 0.05 1.366 H-H!oi 0.31±0.03 

RXJ131 1.5-0551 3.15 ±0.36 0.22!^ ^ 0.692 3.10 ±031 0.20!^ 

RXJ0522.2-3625 3.28!°" 0.14!° : ^ 0.703 3.88!°" °- 31 -fii6 

RXJ2359.5-3211 3.35!° : ^ 0.41!°,^ 0.708 3.26!° : | 0.31!° : |* 

CLJ0030+2618 4 -69!^ 0M -t m °' 836 5 - 01 -032 0.43 ± 0.09 

WARP J2 146.0+0423 3.64!°,^ 0.47!° : fs 0.717 3.40!™ 0.50!°^ 

RXJ0018.8+1602 3.69!° : ^ 0.45!° : ^ 0.718 3.99!°^ 0.60!° : ^ 

MS0015.9+1609 9.48!° : ^ 0.19 ±0.05 1.184 9 - 9 % + nll 0.23 ± 0.05 

WARPJ0848.8+4456 121 -in 0Al+ -afi °- 407 1 - 28 -fll °- 08 -oo8 

CLJ 1354-0221 ZlS -Qm < 0-40 0.542 2.71!°g 0.47!° : ^ 

WARPJ1419.9+0634E 2.86!° : ^ °- 29 !n?s °- 624 2 - 93 to23 °- 32 qw 

MS0451.6-0305 9 - 10 -Q53 0A4 -m L153 9 - 26 -039 °- 18 -qq5 

RXJ0018.3+1618 4 - 65 -»l °- 31 -oi6 °- 807 3 - 87 -Q|l °- 47 -022 

WARPJ1419.3+0638 3 - 48 -o4i < 0.15 0.683 3.93!°™ ao3 -op3 

MS2053.7-0449 3 ' 94 -ol °- 70 -033 °- 725 4 - 44 -058 °- 57 -024 

MACSJ0647.7+7015 9 - 30 !n37 0.27 ± 0.05 1.138 10.52 ± 0.40 0.37 ± 0.04 

CLJ1 120+43 18 5.26!™ 0.34!°!^ 0.837 5.22!°'^ °- 52 -°c,n 

CLJ1334+5031 4 - 19 -04? °- 16 -o!2 °- 733 4 - 73 -Q4i °- 15 -ow 

MACSJ0744.9+3927 8.14!°fj 0.29 ± 0.04 0.995 8.02!°| 0.32 ± 0.03 

WARPJ 1342. 8+4028 3.68!° : " 0.68!°^ 0.652 3.60 ±0.47 0.58!°^ 

CLJ1 103.6+3555 4 - 11 -q| < a4 ° a659 4 23+ °QU a04 -««t 

MS1137.5+6625 6 - 12 -ol °- 19 -SU? °' 813 6 - 60 -«93 °- 37 -on 

CLJ1216.8-1201 3.84!° : !? 0.26!° : ^ 0.631 4.40!° : S 0.09!°' 



MS1054.4-0321 8 - 92 -080 °- 15 -»!r 0969 8 ' 85 -fl! 0.17!"™ 

CLJ0152.7-1357 6.70!™ 0.23!° : |jj 0.829 6.58!™ °- 24 -on 

CLJ1226.9+3332 10 -16!o" 0.21 ±0.08 0.998 11.15 ± 0.58 0.21!°* 

CLJ1429.0+4241 4 - 10 -Q77 < °- 12 °- 606 4 - 75 ^'™ < 010 

XLSSC029 3 - 4 0!o4o 0>38 3na a509 4 - 06 -042 0.56!°^ 

RDCSJ1252-2927 4 -64!° : ^ 0.24!° : |« 0.540 4.84!°|' °- 33 5u 

1WGAJ2235.3-2557 5.16!^ : °° Q-14! ,^ 0-524 5.87^ : l% 0.12!°^ 



To compute the global temperature (T) necessary to esti- 
mate 7-500, we extracted spectra including emission going from 
0.15r5oo to 0.5r5oo in each cluster. The central regions of each 
cluster are therefore excluded from the spectra in order to avoid 
the contamination of a possible cool core. The values of (T) and 
7-500 were evaluated iteratively until a convergence to a stable 
value of the temperature was obtained (Ar < 0.01 keV between 
two successive iterations). Figure [2] shows the distribution of 
< kT > in our sample, ranging from ~ 1 .5 keV to ~ 1 1 keV 
and peaking around ~ 4 keV. 

From the fits, we were also able to determine a global metal- 
licity (Z) in each cluster. In Table |2j we list the best-fit values of 
(T) and (Z), and the value of rsoo computed using the formula 
above. For completeness, in Table [2] we also report the values 
of the temperatures and abundance measured including the core 
emission. 

3.2. Background treatment 

The background of our MOS spectra was modeled (in- 
stead of subtracted) following the procedure developed by 



iLeccardi & Moleridl (1200 8b). These procedure is more accurate 
then a direct background subtraction because it allows us to use 
Cash statistics, which is recommended especia lly for low S/N 
spect ra similar to those in most of our sample dNousek & Shud 
1989), and it enables us to avoid problems caused by the vi- 
gnetting of the background spectra, which need to be extracted 
at larger off-axis angle than the cluster spectra (6 ~ 5') where the 
vignetting of XMM mirrors is not negligible (Fig [3]). Moreover, 
we note that this approach represents the most reliable character- 
ization of the background (and its errors) available with respect 
to any other existing method. The main aspects of this procedure 
are described below, and further details can be found on their pa- 
per. 

We first extracted a spectrum in an annulus located between 
10' and 12' from the center of the field of view for both MOS1 
and MOS2. To estimate the background parameters in this re- 
gion (free of cluster emission), the resulting spectrum was fitted 
in XSPEC (in the 0.7-8 keV band, applying Cash statistics) us- 
ing a model considering thermal emission from the Galaxy halo 
(HALO, XSPEC model: apec), the cosmic X-ray background 
(CXB, XSPEC model: pegpwrlw), a residual from the filter- 
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Fig. 1. Redshift distribution of the XMM-Newton galaxy cluster 
sample. 
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Fig. 2. Average temperature distribution (measured iteratively in 
the 0.15-0.6r5oo range) for the clusters in our XMM-Newton 
sample. 



ing of quiescent soft protons (QSP, XSPEC model: bknpower), 
the cosmic-ray-induced continuum (NXB, XSPEC model: bkn- 
power), and the fluorescence emission lines (XSPEC model: 
gaussiari). The response matrix file (RMF) of the detector was 
convolved with all the background components, while the an- 
cillary response file (ARF) was convolved only with the first 
two components (HALO and CXB). The normalization of the 
QSP component was fixed at the value determined by measur- 
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Fig. 3. Vignetting expected at different off-axis angle for the two 
MOS detectors. The red, green, and blue lines represents the 0.7- 
1 keV, the 1-2 keV, and the 2-8 keV bands, respectively. The 
dashed black lines represents the typical off-axis angle where 
we would extract the background for our sample. It is clear how 
the background spectra could be vignetted up to a 20% factor. 



ing the surface brightness S B m in the 10'- 12' annulus, and com- 
paring it to the surface brightness SB ou , calculated outside the 
field of view in the 6-12 keV energy band. Since soft protons 
are channeled by the telescope mirrors to within the field of 
view and the cosmic -ray-induced background covers the whole 
detector, the ratio R$b = is a good indicator of the in- 

tensity of residual soft protons and was used in the modeling 
(Norm(QSP) = 0.03(1 - R SB )). We also determined the lcr er- 
ror for the background parameters to be used in fitting the clus- 
ter spectra. To fit the cluster spectra, the background parameters 
(and their lcr errors) were then rescaled by the area from which 
the cluster spectra were extra cted. Appropriate correction factors 
dLeccardi & Molend i 2008b, dependent from the off-axis angle) 
were considered for the HALO and CXB components and a vi- 
gnetting factor (corresponding to 1.858 - 0.078r, where r is the 
distance from the center of the annulus) was applied to the QSP 
component. All these rescaled values (and their lcr errors) were 
put in a XSPEC model that had the same background compo- 
nents as those used in the fit of the 10'-12' annulus plus a ther- 
mal mekal model for the emission of the cluster whose temper- 
ature, abundance and normalization were allowed to vary. The 
lcr errors in the parameters were used to fix a range where the 
normalizations of the background components were allowed to 
vary. This model was used to fit all the cluster spectra in our 
sample. Since we used Cash statistics, the spectra were mini- 
mally grouped to avoid spectral bins with zero counts. Moreover, 
a joint MOS1 plus MOS2 fit was used to increase the statistics. 
We were able to do that since there are no calibration problems 
between the two detectors and an independent analysis of the 
spectra from MOS 1 and MOS2 led to consistent results in the 
measures of both kT and Z. 
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Table 3. Comparison of fit results obtained with I Anders & Grevessd (Il989t) (ANGR89) and lAsplund et alJ (120091) (ASPL09) solar 
abundance standards for four of the brightest clusters in our sample. 



Cluster 


Radial bin 

0'500) 


ANGR89 
ytr(keV) Z(Z G ) C-stat 


ASPL09 
ytr(keV) Z(Z Q ) C-stat 


KXJzzz8.0+zl)3 / 


U-U.lj 

0.15-0.4 
0.4-0.9 


8.3/ j^j u.z / Q7 /ui.y 
7.88 + ||f 0.27 + oos 899.4 
5.74+jj^ 0.12 ±0.07 1032.4 


8.58 0.37 010 702.4 
7.89 + |j| 0.38 + °;f 899.2 
5.72+° : :" 0.17 + ° : }' 1032.4 


KXLJ lzUo.z-Uo4o 


U-U. lj 
0.15-0.4 
0.4-0.7 


9.91 _„,, 0.28_ 005 878.2 
11.61™| 0.18 ± 0.07 1120.6 
9.20+]f 0.50^-H 1109.7 


O QC+U.51 r\ OQ+U.OS Q-7Q c 

^■^-047 U '^-Q07 o/o-J 

H-6lt 0.26S H20.6 
9 - 15 -0 82 0.72 ±0.21 1109.1 


RXJ1347.5-1145 


0-0.15 
0.15-0.4 
0.4-0.8 


11.11 ±0.26 0.34+""' 1204.2 

11.55+°$ °- 25 QQ5 1072 - 5 

7.74+<> : || 0.15+™ 1075.3 


11.14 ±0.26 0.48 ±0.05 1205.0 
0.35+°$ 1073.1 
7.57+™ 0.21+° : j 1076.6 


MS0451.6-0305 


0-0.15 
0.15-0.4 
0.4-0.7 


9.47+™ 0.24+™ 631.3 
9.47+° : | 0.18 ±0.07 786.8 
7.29+™ 0.07+°$ 806.2 


9.49+"^ 0.34+°^ 631.2 
9.46+ : | 0.25+°^ 786.7 
7.25+™ 0.10!«^ 806.3 



3.3. Effects of the XMM-Newton PSF on the analysis 

In contrast with Chandra, the EPIC XMM-Newton telescopes 
have a quite large PSF (HEW ~ 15"). This was a potential prob- 
lem to our analysis because the emission from the cluster cores 
(r < 0.15r5oo) in principle could contaminate the spectra just 
outside the core (0.15r5oo < r < 0.4r5oo) and viceversa, if the 
extraction region for the cores is similar in size to the HEW or 
smaller, as happens in about a third of our sample and in almost 
all the clusters at z > 0.7. However, the contribution of the emis- 
sion from the 0.15 - 0.4r5oo annulus to the outer region spectra 
(r > 0.4r5oo) can be assumed to be negligible. If the source is 
uniform over the field, then this effect is not a major concern. 
However, if there are strong gradients in the emission and spec- 
tral parameters over the field (as it is the case of cooling core 
clusters) then the effect might be significant. 

To quantify the effect of the PSF on the determination of 
the spectral parameters measured (and especially on the abun- 
dance measurements), we us ed a procedure developed in XMM- 
ESAS dSnowden et alj |2004) that considers the crosstalk ARFs 
between two contiguous regions to remove the contribution of 
the core from the surrounding region and viceversa. We ap- 
plied it to a few clusters in the sample (RXJ1213.5, CLJ1354, 
MS1054.4, and CLJ1429.0) at different redshifts and core ex- 
traction radii, which shows evidence of a temperature dip at the 
center and abundance gradients. In all of these clusters, the dif- 
ferences in the temperatures measured considering the PSF cor- 
rection and the temperatures measured when not considering the 
correction were in the range between 5% and 15%, which are 
comparable to or smaller than the statistical errors in the mea- 
surements. On the other hand, the differences in the abundance 
measures were always < 0.02 Z , i.e. well below the statistical 
errors. Therefore, we decided not to apply the correction to the 
remainder of the sample. 



3.4. Spectral fitting 

The spectra were analyzed with XSPEC vl2.5.1 (lArnaudl 
fl996h and fitted with a single-temperature mekal model 
( Liedahl et al J 1 1 995b with Galactic absorption (wabs model), 
in which the r atio between the elements was fixed to the so- 
lar value as in lAnders & Grevessel d!989l) . We prefer to report 
metal abundances in these units of solar abundances because 
most of the literature still refers to them. We tested the robustness 
of our results to variations in the assumed emission abundance 
model by fitting the data for four of the brightest clusters in 



our sample (RXJ2228.6+2037, RXCJ1206.2-0848, RXJ1347.5- 
1 145, and MS 45 1.6 -0305) using also the most recent values of 
Asnl und et al.l (l2009h . Table [3] shows the results of the fits ob- 
tained using the two different solar abundance references. The 
fits were of similar quality, as can be seen from the similar values 
of the C-statistic, and the only significant differences were those 
found in the values of the fitted abundances, which we found to 
be consistent with a sim ple scaling by a factor of 1.41 (within 



a 3% margin), where the Asolund et al. (2009) abundances are 



(p_ 

higher than those o f Anders & Grevesse ( 1989). This is a clear 
indication that we have measured the iron abundance in the ICM, 
as expected from the statistics we have in virtually all the cluster 
spectra. In addition, also the fitted temperatures using the two 
abundance models were all consistent. 

All the spectral fits were performed over the energy range 
0.7 - 8.0 keV for both of the two MOS detectors. We fitted si- 
multaneously MOS 1 and MOS2 spectra to increase the quality 
of the statistics. 

The free parameters in our spectral fits are temperature, 
abundance, and normalization. Local absorption is fixed to the 
Galactic neutral hydrogen co lumn density, as obtained from ra- 
dio data (Kalb erla et alJ2 005). and the redshift to the value mea- 
sured from optical spectroscopy (z in Table [TJ- We used Cash 
statistics applied to the so urce plus backgro unqj, which is prefer- 
able for low S/N spectra (Nousek & Shud [T989) and requires a 
minimal grouping of the spectra (at least one count per spec- 
tral bin). We estimated a goodness-of-fit for the best-fit mod- 
els using the XSPEC command goodness, which simulates a 
user-defined number of spectra based on the best-fit model it- 
self and calculates the percentage of these simulations in which 
the fit statistic is smaller than that for the data. For each spec- 
trum in the sample, we simulated 10000 spectra from the best-fit 
model, whose C-statistic values were compared with the value 
obtained from the original spectrum. For the spatially resolved 
bins (r = - 0.15r 500 , r = 0.15 - 0.4r 500 and r > 0.4r 500 ), 
we found that the percentage of simulations with a smaller value 
of the C-statistic than that for the data is always less than 40%, 
an indication that the one-temperature model is always a good 
fit to the spatially resolved spectra. For the spectra extracted in- 
cluding emission from the whole cluster (r = - 0.6r500 and 
r = 0.15 — 0.6r500), the percentage in a few cases is found to be 
larger than 40%, especially in the spectra with more counts. This 
is fully expected because the potential well of galaxy clusters is 



2 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/ 
XSappendixStatistics.html 
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Fig. 4. Left: Abundance measured by EPIC XMM-Newton as a function of the redshift for the clusters in our sample including the 
emission from the whole cluster at r < 0.6rsoo (top) and excising the core (bottom). The error bars in Z are at lcr, while shaded areas 
show the weighted mean of the abundance with its error (blue) and rms dispersion (cyan) in three redshift bins. Right: Abundance 
measured by EPIC XMM-Newton as a function of the redshift for the clusters in our sample in different spatial bins (from top to 
bottom: 0-0.15r5oo,0.15-0.4r5oo, and > 0.4r5oo), with lcr error bars indicated. 



very likely to harbor a multi-temperature gas. In those cases, the 
one-temperature model is not a good fit to the data, but this is not 
a major problem because we simply wish to measure an average 
temperature and abundance for the integrated emission. 

4. Results 

The main results we attained from the spectral analysis are pre- 
sented in this Section. In particular, we focus on the redshift 
evolution of the metal abundance considering both the integrated 
and spatially resolved emission from the clusters, and on the cor- 
relation of the temperature with both redshift and abundance. 

4. 1. Integrated abundance evolution 

Although the main focus of this paper has been to spatially 
resolve the evolution of the abundance with redshift, a com- 
plementary measurement is an integrated measure of the abun- 
dance, which is useful when comparing with previous works in 
the literature. 

Figure [4] (left) shows how the measured abundance in our 
cluster sample evolves with redshift. To perform a statistical ver- 
ification of this evolution, we computed a Spearman's rank p 
for the distributions plotted in the graph considering both emis- 
sion from the whole cluster (r = — 0.6rsoo) and after excis- 
ing the cores (r = 0.15 - 0.6r5oo). We performed 10000 Monte 
Carlo (MC) simulations to take into account the distribution of 
the errors in Z in the determination of the Spearman's rank p. 
In the distribution considering emission from the whole clus- 
ter, we observed a weak negative correlation between redshift 
and abundance with a mean Spearman's rank p = -0.12 + 0.10 
with 39 degrees of freedom (d.o.f.), corresponding to a mean 
probability of no correlation p = 0.48 ± 0.32. The negative cor- 



relation was found to be even weaker after we had excised the 
core, as could be deduced from the computed mean Spearman's 
rank p = -0.05 ±0.11 (again with 39 d.o.f.), which corresponds 
to a probability of no correlation p = 0.61 + 0.31. We also fit- 
ted the data points in Figure [4] (left) with a power-law such as 
Z oc (1 + zY y ', where Z is the abundance measured and z is the 
cluster redshift. The value of y (reported in Table [4]i indicates 
a weak negative evolution with a significance lower than 2cr in 
both cases. The most significant deviation from y — is ob- 
served in the - 0.6rsoo sample, where y = 0.75*912 at a 90% 
confidence level. In the sample where the core is excised, the 
power-law slope is consistent with zero (y = 0.61*21?). 

4.2. Spatially resolved abundance evolution 

As stated before, the main goal of our paper has been to study 
the evolution of the metal abundance in galaxy clusters through 
cosmic time, trying to spatially resolve the cluster emission into 
three radial bins. Figure|4](right) shows how the abundance mea- 
sured at different radii in the cluster sample evolves with red- 
shift. 

Similarly to the steps taken in § 14.11 we used two differ- 
ent statistical methods to quantify any presence of evolution: the 
calculation of a mean Spearman's rank p (averaged performing 
10000 MC simulations) and a power-law fit to the data points. 

In the r < 0.15r5oo radial bin, we do not observe any ev- 
idence of abundance evolution with a mean Spearman's rank 
p = 0.03 + 0.14 with 26 degrees of freedom (d.o.f.), correspond- 
ing to a mean probability of no correlation p = 0.62 + 0.28. 
Fitting the data points in this bin with a power-law in the form 
Z oc (1 + z yy, as in §gj] we obtained y = -0.19+°;^ (Table©, 
indicating no significant evolution in agreement with the mean 
Spearman's rank result. 
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Fig. 5. Temperature vs. Abundance measured for the galaxy 
clusters in our sample including (top) or excluding the core (bot- 
tom). The dotted lines represent the lcr errors in kT and Z, while 
shaded areas show the weighted mean in the abundance with its 
error (blue) and rms dispersion (cyan) in three temperature bins. 
Althought the values of Z show a larger spread at kT < 5 keV, 
we do not find any correlation between the two quantities. 



Fig. 6. Redshift vs. temperature measured for the galaxy clusters 
in our sample including (top) or excluding the core (bottom). 
The errors in kT are plotted at lcr, while shaded areas show the 
weighted mean of the temperature with its error (blue) and rms 
dispersion (cyan) in three redshift bins. No clear correlation be- 
tween redshift and temperature is observed. 



In addition, at 0.15rgoo < r < 0.4r5oo, the values are 
very widely dispersed as is clear from the computed mean 
Spearman's rank p = -0.06 + 0.12 with 39 d.o.f., which cor- 
responds to a probability of no correlation p = 0.54 ± 0.31. 
In this case, the value of y from the power-law fit is positive 
(y = 0.44*22?, as shown in Table @J, although it is consistent 
with zero within the lcr errors, confirming that no significant 
abundance evolution was found in this radial bin. 

The correlation is not present at larger radii (r > 0.4r5oo) 
with a mean Spearman's rank p = 0.08 + 0.16 and a mean p = 
0.54 + 0.33. Fitting the data points with a power-law, the value of 
y is positive but also consistent with zero within lcr in this case 
(y = 1.65+j as reported in Table|4|i. 

4.3. Temperature-abundance and temperature-redshift 
correlation 

As stated in § 13.31 the PSF of EPIC could in principle have an 
effect on the determination of the temperature in the central two 
annuli of the clusters (especially at high redshift) of the order 
of 5-15%, if there were temperature gradients. Therefore we de- 
cided to compare the distribution of Z as a function of the tem- 
perature only for the whole cluster emission (or after simply ex- 
cising the core as in § 14. U . 

Figure [5] shows the correlation between the temperature and 
the abundance measured in each cluster including (r = - 
0.6r5oo) and excluding the core (r — 0.15 - 0.6/500). Although at 
first glance, it appears that higher values of abundance are more 
frequent at low temperatures, the correlation between this two 
quantities is statistically very weak (as can be clearly seen from 
the weighted average values of Z derived in three different bins 



of temperature). When we considered the emission from the 
whole cluster, we measured a Spearman's rank p = -0.09 + 0. 1 3 
with 39 d.o.f. (errors computed using 10000 MC simulations 
as in § 14. Q , corresponding to a probability of no correlation 
p = 0.52 + 0.33. The correlation was also found to be weak 
after we had excised the core (p = 0.02 + 0.13 with 39 d.o.f., 
p = 0.56 + 0.30). We also fit the distribution with a power-law of 
the form Z oc kT~ y . The results are consistent with there being no 
correlation between kT and Z for both samples (y = -0.06 +0. 1 1 
including the core, and y = -0.06 + 0.16 excising the core). 
Since, as stated before, we found that the values of Z have a 
greater dispersion at lower temperatures, we divided each sam- 
ple into two subsamples at kT < 5 keV and kT > 5 keV and 
compared them with a Kolmogorov-Smirnov (K-S) test to verify 
whether they come from the same distribution. Both including 
and excising the core, the K-S test cannot reject at a 5% con- 
fidence level the null-hypothesis that the clusters with kT < 5 
keV have the same distribution of Z as the clusters with kT > 5 
keV. We therefore do not find any obvious trend in the abun- 
dance with temperature similar to those seen in other samples 
of galaxy clusters observed with other X-ray missions suc h as 
Chandra (BLS07) and ASCA (e.g. lBaumgartner et alJl2005h . 

No clear correlation is also observed in the redshift ver- 
sus (vs.) temperature distribution (Figure |6), as indicated by a 
Spearman's rank p = 0.14 + 0.11 with 39 d.o.f., (probability of 
no correlation p = 0.39 + 0.25), for the - 0.6rsoo sample, and a 
Spearman's rankp = 0.09 ±0.11 (39 d.o.f., p = 0.51 + 0.26) for 
the 0.15 - 0. 6/500 sample. Figure|6]also shows that we have ana- 
lyzed a population of clusters whose temperatures are uniformly 
distributed with redshift, there being no preferential redshift 
range containing the hottest clusters in the sample (kT > 8 keV). 
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Table 4. Values of y in the power-law fits of the z-Z distribution (Z oc (1 +zT y ) in the samples present in the literature compared with 
those obtained in the present work. The errors in y and the normalization (referred at z = 0.6) are quoted at lcr for one interesting 
parameter. 





Extraction radius 


Nobj 


Redshift range 


r 


Norm (@z = 0.6) 


inis worK < q|4. tD 


- U.6r 50 o 

ft 1 S ft fir,,.,. 

U.l J U.OA500 


on 

jy 
jy 


U.4U/— l.jyi 

U.+U f—l .J7J 


n 7 c +0.48 
06+^ 


U.Zo ± U.U1 
OA -+- O 09 


This work fij 14.21 


0-0.15r 500 
0.15-0.4r 500 
> 0.4r 500 


26 
39 
24 


0.407-0.920 
0.407-1.393 
0.407-1.237 


-°- 19+ a68 
44+°™ 

?:S+?f 


0.33 ± 0.02 
0.23 ± 0.02 
0.13 ±0.03 


BLS07 


0.15-0.3r„„. 


46 


0.405-1.273 




0.28 ±0.01 


MAU08" 


- lr 50 o 
0.15-lr 500 


50 
46 


0.405-1.237 
0.405-1.237 


7 44+8.™ 

q 07+ 1.06 


0.30 ± 0.02 
0.22 ± 0.02 


MAU08" 


- lr 50 o 
0.15 - lr 500 


49 
45 


0.405-1.237 
0.405-1.237 


2.18_ 07 , 
-0 15 +lfe 


0.31 ±0.02 
0.28 ± 0.02 


Anderson et al. (2009) 


- X" c 


23 


0.407-1.237 


09+ m 


0.24 ± 0.02 



Notes. {a) Including CLJ1415. 1+3612 in the fit. 
(W Not including CLJ1415. 1+3612 in the fit. 

(c) Maximum extraction radius X listed in their paper and determined using S/N optimization criteria. 



5. Discussion 

The main goal of our paper has been to study the evolution of 
the metal abundance in galaxy clusters through cosmic time. To 
achieve this aim, we did not consider the emission only from 
the entire cluster (§ 14. Il l but we also tried to spatially resolve 
the cluster emission into three radial bins of - 0.15r5oo, 0.15 - 
0.4r 500 , and>0.4r 50 o(§i3- 

5. 1. Comparison with previous works 

In § I4.U we have shown that considering the integrated emission 
from the whole clusters we could detect only a mild negative 
abundance evolution with redshift at a confidence level lower 
than 2<x (as indicated by the power-law fits of the z-Z distribu- 
tion). 

Notwithstanding, it is interesting to perform a quantitative 
comparison of our results with those obtained previously in 
the literature (e.g. BLS07). The Chandra sam ples of BLS07 
and M AU08 and the XMM-Newton sample of I Anderson et al.l 
(2009) are the most suitable to perform this comparison. In these 
samples, we selected only the clusters at z > 0.4 to be consis- 
tent with the redshift range probed in our work. We have 23, 46, 
50, an d 46 objects in this redshift range for the An derson et al.l 
(l2009l) sample, the BLS07 sample, the MAU08 sample, and the 
MAU08 sample with the core excised (r < 0.15r5oo), respec- 
tively. We performed the same power-law fit to the z-Z relation 
that we used in § l4.1l for our data (in the form Z oc (1 +z)~ y ). The 
comparison of the fits with those obtained for our sample are 
shown in Table [4] and Figure [7] where the confidence contours 
on y and normalization (referred at z=0.6) are shown for all the 
datasets analyzed. From the values of y and normalization shown 
in the table and the confidence contours shown in Figure [7] it is 
clear how our results for the integrated measure of the abundance 
in the - 0.6r5oo sample (ij l4.ll ) are fully consistent with BLS07 
sample. Both these samples do not display any significant signs 
of negative evolution (sig nificance lower than 2c ). Table [4] also 
shows that the sample of lAnderson et al.l d2009l) has a y ~ 0, 
while the only samples showing an evolution of abundance with 
z at a confidence level higher than 2cr are the MAU08 Chandra 
samples. In particular, it is striking how the values of y from 
the MAU08 sample are significantly higher than in all the other 
samples (especially not considering the core emission). We note 



that these values are probably affected by the presence in their 
sample of a single cluster at high redshift with a very low up- 
per limit measurement on Z (CL J 1415. 1+3612 at z - 1.030 with 
Z < 0.14 and Z < 0.04, including and excising the core, re- 
spectively). Therefore we performed the power-law fits of their 
sample excluding this particular cluster to quantify its influence 
on the z vs. Z trends. While in the sample including the core, 
the results of the evolution strength are quite robust (y changes 
from 2.44+° ™ to 2.18+°^), removing CLJ1415. 1+3612 from the 
0.15 - lrsoo sample has a dramatic effect on the measured abun- 
dance evolution. The value of y indeed changes from 3.87+igj, 
corresponding to a strong negative evolution, to -0.15*H?, cor- 
responding to no significant evolution. This is also evident from 
Figure [7] where the confidence contours in the y-normalization 
parameter space referring to this sample are consistent with no 
evolution and with all the other datasets considered in this sec- 
tion. That a single deviant object could have such a strong effect 
on the abundance evolution significance is clearly a caveat about 
being careful in interpreting results based on such small sam- 
ples. In any case, from the values reported in Table [4] and the 
projection of the contours in Figure|7| we can constrain the nor- 
malization (at z = 0.6) to Z as 0.28 andZ « 0.24inour0-0.6r5oo 
sample and in our 0. 15 - 0.6r5oo sample, respectively. These val- 
ues are consistent within lcr with the values obtained from all the 
other datasets in the literature considered in the present section. 
We note that the lack of evolution in almost all the samples ana- 
lyzed in this section could be due, at least in part , to the red shift 
range probed in this paper. BLS07, MAU08 and And erson et al.l 
(2009) did indeed detect a significant evolution in abundance, 
only between the low redshift clusters (0.1 i; z 5 0.3) and the 
clusters at z > 0.5. 

As a consistency check, we also applied a Bay esian MC 
metho d to fit the datasets described above (see e.g. lAndreonl 
2010), finding consistent results with the power-law minimum 
least squares fit for all the samples. 

5.2. Spatially resolved abundance: weighted mean and 
radial dependence ofz 

In § 14.21 we resolved spatially the emission from our clusters into 
three radial bins. However, no evolution in abundance could be 
observed in each of the radial bins considered. In this case, the 
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Table 5. Mean weighted abundance in the three radial bins considered computed at three different redshift bins, with their errors 
and r.m.s. dispersions indicated. 



Redshift 


<Z(0 - 0.15r 50 o)> ± err ± rms 


<Z(0.15 - 0.4r 50 o)> ± err ± rms 


(Z(> 0.4r 50 o)> ± err ± rms 


0.4-0.5 
0.5-0.7 
0.7-1.4 


0.32 ±0.02 ±0.19 
0.37 ± 0.03 ± 0.29 
0.28 ± 0.09 ± 0.20 


0.25 ±0.03 ±0.17 
0.27 ± 0.03 ± 0.20 
0.19 ±0.05 ±0.08 


0.17 ±0.04 ±0.16 
0.15 ±0.04 ±0.23 
0.09 ± 0.08 ± 0.44 




Fig. 7. Confidence contours (computed at lcr and 2<x for two in- 
teresting parameters) on y and normalization (referred to z=0.6) 
for all the datasets analyzed in § 15.11 Blue and cyan contours 
refer to our - 0.6r5oo sample and our 0.15 - 0.6r5oo sample, 
respectively. The contours for MAU08 samples (computed ex- 
cluding CLJ1415. 1+3612, see text and Table© at r = - lr 500 
and r = 0.15 - lrsoo are plotted in red and pink, respectively. 
Green contours refer t o the BLS07 sample, w hile orange con- 
tours are plotted for the I Anderson et all d2009l) dataset. 



deviations from no evolution are even less significant than in the 
integrated measures we have obtained in this paper (§|4Tj and 
in the results of BLS07 and MAU08. The most likely reasons 
for this discrepancy are the larger errors in the measure of Z and 
the smaller sample size of the r = - 0.15rsoo and the r > 
0.4r5oo samples (26 and 24 clusters, respectively). Both of these 
reasons are clearly caused by the smaller number of counts we 
have when dividing the emission into three radial bins instead of 
considering a single spatial bin for each cluster. 

Although the low quality of our statistics prevent us from 
any quantitative claim about the abundance evolution, it is in- 
teresting to consider a weighted mean of the abundance com- 
puted at different redshift bins to perform a qualitative compar- 
ison with the literature (in particular wi th MAU08) and the the- 
oretical predictions of the Ettori (2005) model. The size of our 
XMM-Newton sample did not allow us to consider more than 
three redshift bins (0.4 < z < 0.5, 0.5 < z < 0.7, and z > 0.7) 
when computing the weighted mean and using as weights the 
statistical errors in the measure of the abundance. Figure[8]shows 
the mean weighted abundance in the three radial bins considered 
computed at the three different redshift bins listed above. The 
mean weighted abundances are also shown in Table|5]with their 
errors and r.m.s. dispersions. We have also plotted the weighted 
means of MAU08 for their 116 Chandra clusters (obtained by 
performing a stacking spec tral analysis a nd excising the cores) 
and the predictions of the lEttoril d2005h model for the abun- 



0.4 



0.3 



, 0.3 



0.1 




Redshift 

Fig. 8. Mean weighted abundance in our XMM-Newton sample 
as a function of the redshift for the clusters in our sample in 
different spatial bins: 0-0.15r5oo plotted in red, 0.15-0.4r5oo in 
green, and > 0.4rsoo in blue. The error bars represent the error in 
the mean and are listed in Table The cyan points represent the 
weighted abundance measured by MAU08 for a sample of 116 
clust ers. The ab undance values expected at different redshifts in 
the lEttoril (l2005h model are plotted as a solid black line. 



dance evolution in galaxy clusters. As expected, the abundance 
weighted-mean decreases towards larger radii (consistently with 
the power-law normalizations reported in Table©). However, the 
z-Z relation is consistent with no evolution at every radius. As 
stated before, the small number of counts in each radial bin (and 
the associated large errors in Z) for most of the clusters in the 
sample does not allow us to measure with sufficient accuracy the 
mean abundance at different redshift, and prevents any statisti- 
cally significant claim about its evolution with cosmic time in 
the three different cluster regions considered (i.e. the core, the 
core surroundings, and the outskirts), in agreement to what we 
found in § 14.21 From Figure [8] we can see that, although there is 
no supporting statistical evidence of any kind of evolution, the 
trend in the 0. 15-0.4r5oo radial bin generally complements nicely 
the measures of MAU08, obtained after excising the cool core 
(0. 15r5oo < r < rspp), and broadly agrees with the predict i ons of 
the lEttoril (120051) model. In particular, the model of lEttoril (120051) 
predicts a slope between z = 0.4 and z - 1.4 of y = 1.7, which 
is only ~ 1.5<x above the constraints we have in our 0.15-0.4r5oo 
dataset. 

The dependence of the abundance Z on both the radius r and 
the redshift z could also be tested by fitting Z, r, and z from the 
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89 data points in the three radial bins together (the points plotted 
in the right panel of FigureH]) with a function of two variables in 
the form 



Z(r,z) = Zo 



1+1^ 



2\ 



ill 
1.6 



(2) 



In this formula, the radius of reference r in each radial bin was 
defined as the effective radius R weighted by the emissivity dis- 
tribution within a given radial bin defined by 



2 f fi ° 
S b (R h R )= , S b rdr = S 



(3) 



where Sb is the X-ray surface brightness, and and R Q are the 
inner and outer radius of the radial bin. Under the assumption 
that S b x T~ a , we obtain 



R 



2 (R 2 a -R*- a )) 
(2-a) X (R1-R*) 



- 1 /a 



(4) 



In our estimates, we adopt a = 0.7, 2.2, and 3 in the - 0.15, 
0.15 - 0.4, and > 0.4r5oo radial bin, respectively, in accor- 
dance with the expected radial behaviour of the X-ray surfac e 
brightness in massive galaxy clusters (see e.g. lEttori et al.l2 004). 
Fixing ro = 0.15rgoo, we obtained a good fit with a% 2 ~ 107.6 
with 86 d.o.f. ix^-ed ~ 1-25). The best-fit values of the three 
free parameters are Zq = 0.36 + 0.03, a = 0.32 + 0.07, and 
y = 0.25 + 0.57. These values represent a significant negative 
trend of Z with the radius and again no significant evolution 
with redshift. That the evolution with redshift is even smaller 
when the radial dependence is taken into account suggests that 
part of the claimed evidence of a negative evolution of the inte- 
grated metallicity can be ascribed to the decrease with radius in 
the measured abundance, for instance when the metal content is 
estimated in relatively smaller regions in the nearby systems. 

Although the constraints on the abundance evolution with 
the redshift are far from being stringent, this is the first time 
that the abundance measured in the ICM of galaxy clusters has 
been parametrized as a function of both the distance from the 
center, and the redshift. This new result should be directly com- 
pared with th e models of diffusion of metals ("e.g. | Ettprill2005t 
lLoewensteinll2006l: ICalura et alJkOOTt ICora et al.ll2008l) and the 



simulations of chemical enrichment (e.glTbrnatore et alj 12 004 ; 
Rebusco et al. 2005; Tornatore et al. 2007; Schindler & Diaferio 
2008; lFabian et alj|2010h in the ICM. 



6. Conclusions 

We have presented a spatially resolved XMM-Newton analysis 
of the X-ray spectra of 39 clusters of galaxies at 0.4 < z < 1.4, 
covering a temperature range of 1 .5 % kT $ 1 1 ke V. The main 
goal of this paper has been to study how the abundance evolves 
with redshift not only by analyzing a single emission measure 
performed on the w hole cluster (as in e.g. BLS07; MAU08; 
lAnderson et al.ll2009h . but also by spatially resolving the cluster 
emission in three different regions: the core region (correspond- 
ing to < r < 0.15r5oo), the region immediately surrounding 
the core (0.15rsoo < r < 0.4r5oo), and the outskirts of the clus- 
ter (r > 0.4r5oo). This represents the first time that a spatially 
resolved analysis of the metal abundance evolution through cos- 
mic time has been attempted. 

The main results of our analysis can be summarized as fol- 
lows: 



- Across the studied redshift range, we have not found any sta- 
tistically significant abundance evolution with redshift con- 
sidering the integrated emission from the whole cluster both 
at r < 0.6r5oo and at r — 0.15 - 0.6rsoo. The most signifi- 
cant deviation from no evolution was found in the - 0.6r5oo 
sample, where a power-law fit to the z — Z distribution gives 
a y = 0.75+°^ at a 90% confidence level. 

- Dividing the emission into three radial bins, we could not 
find any significant evidence of abundance evolution (al- 
ways below lcr). The main reasons for this are probably the 
smaller number of counts (and the associated larger errors in 
Z) and smaller statistical samples for two of the three radial 
bins. 

- We have not found any significant correlation between 
the abundance and temperature similar to that seen in 
other samples of galaxy clusters observed with other X- 
ray missions such as C handra (BLS07) and ASCA (e.g. 
Baumgartner et al Eool) . Using different statistical methods 
(Spearman's rank, power-law fit, K-S test), the significance 
of a kT - Z relation is always less than lcr. 

- We compared our results for the single-emission-measure 
abundances with previous works in the literature. Our re- 
sults are fully consis tent with the BLS07, MAU08, and 
lAnderson et al.l d2009l) samples at z > 0.4, as indicated by a 
power-law fit to the z-Z distributions. The differences in the 
power-law slope y and in the normalizations at z = 0.6 are 
always smaller than 2<x, with the largest discrepancy found 
when comparing with the MAU08 - lrsoo sample. 

- We found that the presence or not of a single cluster 
(CLJ1415. 1+3612) in the MAU08 sample could dramati- 
cally change the abundance evolution measured at z > 0.4 
from a power-law fit, such that in the 0.15 - lrsoo sample, 
y varies from 3,87*1x5 (indicating strong negative evolution) 
to -0.15^ 23 (no significant evolution). The impact that a 
single deviant object could have on the measured abundance 
evolution is a caveat about being careful in interpreting re- 
sults based on such small samples. 

- We computed error- weighted means of the spatially resolved 
abundances in three redshift bins: 0.4 < z < 0.5, 0.5 < z < 
0.7, andz > 0.7. The abundance weighted-mean is consistent 
with being constant with redshift at every radius, given the 
large lcr error bars. A qualitative comparison w ith MAU08 
and the predictions of the model of Ettorjl (120051) shows how 
the trend in the error-weighted mean abundance with red- 
shift in the 0.15-0.4r5oo radial bin complements nicely the 
measures of MAU08 (at 0.15 r^Q < r < r 5oo), and broadly 
agrees with the predictions of Ettori (2005) model. 

- We quantified the dependence of the metal abundance on 
both the radius r and the redshift z by fitting the data points 
available from the spatially resolved analysis with a function 
of two variables Z(r,z) = Z (l + (r/0.15r 50 o) 2 )^(z/0.6)- y , 
with best-fit model values of Z = 0.36 + 0.03, a = 0.32 ± 
0.07, and y = 0.25 + 0.57. This relation depicts a signifi- 
cant negative trend of Z with the radius and, although no sig- 
nificant evolution with the redshift is detected, it represents 
the first time that the abundance measured in ICM has been 
parametrized as a function of both r and z in a way that could 
be direc tly compared with the models of di ffusion of met- 
als (e.g. lLoewen stein 2006t ICalura et al . 2007) a nd the sim- 
ulations of chemical en richment (e.g ISchindler & Diaferiol 
120081: iFabian etaDl2010h in the ICM. 

In general, the low quality of the statistics (and large errors 
in Z) associated with most of the clusters in the sample has pre- 



11 



A. Baldi et al.: An XMM-Newton study of metal abundance in distant galaxy clusters 



vented us fron drawing any statistically significant conclusion 
about the different evolutionary path the abundance may have 
traversed through cosmic time in the three different cluster re- 
gions considered. 

We stress that the results presented in this paper (from both 
our new analysis of the XMM-Newton exposures available at 
z > 0.4 and those already presented in the literature) represent 
the most robust limits that could be derived from the existing 
XMM-Newton and Chandra archival data on the evolution of 
the metal content of the ICM in massive halos. Larger samples of 
high redshift clusters observed in the X-rays and deeper observa- 
tions of the clusters already present in both Chandra and XMM- 
Newton archive would be crucial to place statistically significant 
constraints on the evolution of metal abundances in galaxy clus- 
ters and to provide a robust modelling of the physical processes 
responsible for the enrichment of the cluster plasma during its 
assembly over cosmic time. 
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